Factorizing numbers with classical interference: several implementations in optics 
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Truncated Fourier, Gauss, Kummer and exponential sums can be used to lactorize numbers: for 
a factor these sums equal unity in absolute value, whereas they nearly vanish for any other number. 
We show how this factorization algorithm can emerge from superpositions of classical light waves 
and we present a number of simple implementations in optics. 
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q '■ I. INTRODUCTION 

(N : 

' Factorization of numbers into their prime factors is a hard non-polynomial problem for classical computers. It was 
Shor [l[ who proposed a quantum algorithm which can solve the problem of factorization of numbers on a quantum 
■ computer with a tremendous speedup as compared to a classical computer. A practical demonstration of Shor's 
algorithm has been carried out by factorizing the integer 15 Q, using nuclear magnetic resonance. However, quantum 
computers capable of implementing Shor's algorithm for larger numbers have not been developed yet. 

Several approaches to factorize numbers based on interference of multiple quantum paths have been proposed 



beveral approaches tc 

i— —i . [1, 0, H, H, 0, H, H, [HI LlK ■ Those schemes do not use quantum entanglement and do not capitalize on quantum 
' parallelism. As a consequence, these schemes scale exponentially with the number of digits of the factorized number. 
This is in contrast to Shor's algorithm which requires only a polynomial number of operations. Nevertheless, if the 
interference is implemented in a suitable way in a system, which does the factorization, then one can benefit because 
nature plays the role of a computer. 

As was pointed out by Jones [ljj the proposed techniques for factorization based on Gauss sums [H, 0, 9 unfortu- 
Q-i nately do not provide useful methods to factorize numbers, because a precalculation of the factors is needed for the 
experiment. In spite of that, Gauss sums would be useful if it is possible to avoid explicit precalculation stages of the 
algorithms. 

Physical systems that can implement the Gauss sums must be described by complex numbers. In the present paper 
■ we investigate how truncated Fourier sum and its generalizations, like truncated Gauss, Kummer and exponential 
| sums, could emerge from superposition of several oscillations. Those sums can be used successfully to factorize 
C^) numbers. Due to the wide use of interferences and beats in optics, we shall keep our consideration close to optics and 
' our examples are in wave optics too. 

However, the proposed implementation can be extended to virtually any physical system where superposition 
among several different oscillations appear, from the mechanical pendulum with several degrees of freedom, through 
OO , the atomic and solid state systems and their analogs in quantum mechanics. 
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II. TRUNCATED FOURIER, GAUSS, KUMMER AND EXPONENTIAL SUMS 



■ In order to find the factors of a given number N we use the following truncated sum: 

A M) (o = ^E ex p(- w f)' « 

m=l ^ ' 

where k is an integer and M is the number of terms in the sum. The argument I scans through all integers between 1 
and \fN for possible factors. The capability of the sum of Eq. (flj to factor numbers originate from the fact that for 
an integer factor q of N with TV = ql, all phases in A^\l) are integer multiples of 2-k. Consequently, the terms add 

up constructively and yield A^\l) = 1 . When I is not a factor, the phases oscillate rapidly with m, and A^\l) 
takes on small values. In this interference pattern, larger truncation parameter M leads to better convergency. In 
principle, already the first several terms of the sum are sufficient to discriminate factors from non-factors. Depending 
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on the coefficient k in Eq.(p} we distinguish several important cases: 

Fourier sum for k = 10 
Gauss sum for k = 2 @, 0, 0, 
Kummcr sum for k = 3 
exponential sum for fc = m 
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(2a) 
(2b) 
(2c) 
(2d) 



The use of quadratic phases to factor numbers (Gauss sum) has the advantage of fewer terms needed in the sum 
to distinguish factors from non- factors compared to the linear phase (Fourier sum), which is because of high quasi- 
randomness for the quadratic phase [Io| . In the very same way the Kummer sum, and sums with nonlinear phases of 
higher order, has an advantage compared to the Gauss sum [10l ] . 

Now we consider a system with M different oscillation modes, with frequencies uj m , phases tp m and amplitudes Eo m 



E m (*) = E 0m exp (iuj m t + i(p m ) . 
Using the superposition principle we can write the resulting oscillation as the sum of all oscillations: 
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E(t)=^2 E rn (t) = ^2 Eo m cxp (iu m t + i(p m ) . 



(3) 



(4) 



m— 1 



m— 1 



In the sum of Eq. Q we can vary the parameters E 0m , uj m , (p m and the time t. In the next several sections we will 
show how truncated Fourier, Gauss, Kummer and exponential sums could emerge when we fix three of the parameters 
for all oscillations, while changing the fourth parameter. 
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III. FACTORIZATION USING DIFFERENCES IN TIME DELAY (INTERFEROMETRY) 

A. Mach-Zehnder interferometer 

First we consider the case when the parameters Eo m , uj rn , tp m are equal for all oscillations in Eq. (|3|) 

(5) 
(6) 
(7) 

thus the only parameter that is left not fixed in Eq. ([3]) is the time t. This can be easily realized in optics by 
interferometry, where the individual oscillations describe the electric field for the different arms of the interferometer 
as shown in Fig[TJ 

From Fig. Q]and Eq. (0| we see that we have the following sums of electric fields in the detector 

M 

E = EoY, ex P , ( 8 ) 

m— 1 

where <j) m is the phase accumulated in the m arm of the interferometer due to the difference in travel time through 
each arm. 

Suppose that each arm of the interferometer is with length L, the wave length of the light that we use in vacuum is 
A and the corresponding frequency is u>, let the index of refraction in each arm of the interferometer is different and is 
denoted as n m . Then the phase <p m for the beam that travels through the m-th arm of the interferometer is given as 

<t> m = tmU = — u, (9) 

Cm 

here t m is the time that light travel in the m-th arm of the interferometer to pass length L and c m is the speed of 
light in that arm. The refraction index in arm m of the interferometer is 

n m = — , (10) 

Cm 
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thus 



then the electric filed in the detector is 



L 

uj—n r , 
c 



M ( L \ 
E = CX P ( 27ri— n m J 

m=l ^ ' 
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-n m L, 



Now if the index of reflection in the m arm of the interferometer is 

n m — a + bm k , 

then 



E = E a cxp ( ) cxp ( 2Trim k ^- 



.bL 



The detector registers the intensity, 



\E\ 2 = 



A I 



exp 2mm — 



m— 1 



(11) 



(12) 



(13) 



(14) 



(15) 



Various k gives us a different type of truncated sum (see Eq. ©). The number that we want to factorize is bL, 
the trial factors are A. Each time when the trial factor A is a factor of bL we will observe a maximum signal in the 
detector. The number of the terms in the sum can be controlled by doubling the elements in the interferometer Fig. 
[TJ The numbers that could be factorized in this way are of order L/A ~ ln n"L„ = 10 6 
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B. Pulse train 



Now we consider a train of pulses, where the delay of the m pulse compared to first pulse is given as 

Ua = m k T, (16) 

here m takes the values in = 1, 2, 3...M, while r can be set as a unit of time. 

We consider the case when all pulses have equal amplitudes Eo m = Eq and equal frequencies uj m = u. Then the 
electric field for the m pulse is given by Eq. and reads 

E m = E exp (iut m + iip m ) = E exp (iujm k T + iip m ) (17) 

Let us make a different path way for every pulse in such a way that all pulses hit the same detector at the same time, 
this is equivalent to make ip m = at the place where all pulses collides. Then the intensity that the detector registers 
is a result from the superposition among all electric fields, e.g the sum from Eq. ([4]): 
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I~\E\ 2 = 






Eq exp (2iirm k vt) 






m— 1 




ra—l 





where v = tu/ (2tt). If one chooses the frequency v as the number that we want to factorize (N) and 1/t as a trail 
factor (I), then Eq. fTg)) reduces to the sum from Eq. ([TJ. 



IV. FACTORIZATION USING DIFFERENCES IN FREQUENCIES (BEATS) 

If we now consider a system that exhibits several oscillations with the same amplitude Eq and the same initial 
phases ((p m = 0), but with different frequencies uj m , then the individual oscillations ([3]) arc described by 

E m (t) = E Q cxp (iu m t) , (19) 
uj m = m k uj , (20) 
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the resulting oscillation (0} is 

N 

£ 0) = ]T E ° ex P (-2nim k v t) , (21) 

m— 1 

where uq = ujq/ (2tt). We will observe beats when tvo is a integer which could be used to find the factors of the 
number v§. One physical realization of the above idea could be a light with several high-harmonic generated frequencies 
[TBI [l6| , chosen in the way that they present for example the odd terms in the Fourier sum: 

wo,3wo,5w ,7wo---- (22) 

then in the detector the time of the detection play the role of the test factors and whenever there is a beat we observe 
a maximum of the signal, thus this time is a real factor. 



V. FACTORIZATION USING FARADAY EFFECT 



The last parameter that we can vary in Eq. (J3j is the amplitude of the individual oscillation. For example if we 
work with laser light we can use the different polarization orientations of the electric field. The electric field is a vector 
in the polarization plane, which can be described by complex electrical field. 

Let us consider the case when we have a linearly polarized light pulse, which is split in several parts and each part 
passes different pathways through Faraday cells as shown in Fig. O 

Applying different Faraday rotation angles ip m on each pathway and collecting all of the light at the same place (at 
the detector) the resulting electric field is the superposition: 

E M 

E =JfJ2 ex P ( i( Pm) > (23) 

m— 1 

where Eq is the electrical field amplitude of the initial beam. The relation between the angle of polarization rotation 
due to the Faraday effect (p m and the magnetic field B m in a diamagnetic material [l7| is 



<p m = 2wbLB m , (24) 
where L is the length of each pathway and 2nb is the Vcrdct constant for the material ■ For the amplitude of the 



resulted electric field we have: 



E 

E = -j ex P (2nibLB„ 



M 



M 



(25) 



If we now have a magnetic field B m for the m-th Faraday cell, which is given as: 



B m = B m\ 



then the intensity in the detector is 



M 



— Y exp (2Trim k bLB ) 



(26) 



(27) 



Here the number that we want to factorize is bL, the trial factors are 1/Bq. 



VI. CONCLUSIONS 



We have shown how the factorization algorithm based on truncated Fourier, Gauss, Kummer or exponential sums 
emerges naturally from superpositions of classical light waves. We have proposed a number of simple implementations 
in optics. These implementations can be extended to virtually any physical system where superpositions of several 
different oscillations appear. 

The factorization algorithms discussed in this paper are classical algorithms and thus their complexity scales ex- 
ponentially with the number of digits. If an extension of this algorithm exists in entangled quantum systems, then a 
quantum computing parallelism would be involved with an exponential speedup of factorization. The present solutions 
therefore could be the first step to an alternative quantum factorization algorithm to the famous Shor algorithm. 
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FIG. 1: Four arms Mach-Zehnder interferometer that can be used to factorize numbers by using Fourier, Gauss, Kummer 
or exponential sums. The four arms correspond to four terms in the sum. Repeating the procedure for doubling the arm in 
principle one can increase the terms in the sum as mach as needed to distinguish factors from nonfactors. 
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FIG. 2: Two different path ways of the light that pass through Faraday cells and are used to factorize numbers. 



